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Abstract 

One examines the interaction and possible resonances between supernova 
neutrinos and electron plasma waves. The neutrino phase space distribution 
and its boundary regions are analyzed in detail. It is shown that the boundary 
regions are too wide to produce non- linear resonant effects. The growth or 
damping rates induced by neutrinos are always proportional to the neutrino 
flux and C7|. 
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1 Introduction 



The study of the interactions between neutrinos and plasma waves has received a 
great deal of attention. Early works of Bingham et al. |T], Q claimed that an intense 
neutrino flux passing through a plasma is capable of producing unstable modes of 
electron density waves causing a significant transfer of energy from neutrinos to the 
mantle of a Supernova. If true, this could constitute a sort of realization of the 
Wilson explosion mechanism. However, their description of the weak interactions in 
the neutrino-electron system was not satisfactory. 

More recently || ||, the full Standard Model Quantum Field Theory of elec- 
troweak interactions was applied to establish the dynamics of the excitations of 
the electromagnetic field and electron and neutrino current density distributions 
(Eqs. (27)-(32) of |4]]). From that we derived the modification on the dispersion 
relation of electron density waves due to a neutrino flow. Our analysis of the condi- 
tions of neutrino emission in Supernovae led us to the conclusion that they do not 
satisfy the necessary requirements to generate the unstable waves and growing rates 
predicted by Bingham et al. |], §. Subsequent works ||, [J gave some support to 
this conclusion. However, the controversy seems to persist on this specific issue H. 
We report here a detailed analysis concerning the neutrino phase space distribution 
and implications for plasma waves. 



2 Neutrino induced instabilities 

The dispersion relation of electron density waves, also called plasmons, is modified 
by weak interactions when a neutrino flow passes through the plasma. Let uj p i(k) 
designate the wave frequency as a function of the wavevector k in absence of neutri- 
nos and in the plasma collisionless limit. We assume the background medium to be 
static and spatially homogeneous within the time and length scales characteristic of 
the plasma waves of interest. The momentum distribution is assumed isotropic for 
electrons but not for neutrinos, as they have an almost unique direction at far dis- 
tances from the Supernova core. The electron plasma is considered non-relativistic. 
In these conditions a stream of neutrinos (anti-neutrinos) modifies the dispersion 
relation as follows 0, |J (h — c — 1) 

u 2 -u 2 pl = TAu] 2 pl , (1) 

T = 2G14^ } (2) 
m e E v 

where c' v = 1/2 + 2sin 2 #w ~ 0.96 for v e (z? e ), Gp is the Fermi constant, n e , n u are 
the electron and neutrino number densities, respectively, and 

u 2 \ 2 k 2 E u f f u k 2 -(k- Vl/ ) 



A ' ' k? ) uoln u j d Vv E v (lo - k.v„)2 - (Co* - k2)2/4£2 • ''•>; 
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is a dimensionless quantity. /„ is the neutrino or anti-neutrino distribution function 
in momentum space, J f u d 3 p u = n u , and E v is a typical neutrino energy. Finally, 
Up = Airan e /m e . 

The expression A admits in general a classical approximation because the fre- 
quency and wavenumber are much smaller than the single particle energy E u . In 
fact uj is around the magnitude of the plasma frequency, lj p , and k is limited above 
by the Debye wavenumber, = u p ^m e /T e at a temperature T e ||. In the classic 
limit A is approximated as 



or, after integrating by parts, 



V k 2 J u l p n v J u}-k-v v 



(5) 



an expression directly related to classic kinetic theory || |]. 

The frequency shift due to weak interactions is in general extremely small. The 
factor G\n e n v jm e E v is about 3 x 1CT 28 for electron and neutrino densities as high 
as n e = 10 29 cm -3 and n u = L u / Airr 2 E u = 1.8 x 10 30 cm -3 at radius r = 300 km 
for a neutrino energy luminosity L v = 10 52 erg/s and E v = lOMeV. The claim has 
been |TJ, ^ |7|] that the shift on the imaginary part of the frequency, 7 = Im{u — ui p i}, 
is not suppressed by Gp but rather by a smaller power of Gf for certain wave modes 
that are resonantly enhanced by powers of (to — k-v^) -1 ~ 7 _1 3> cu^ 1 . That would 
be the case if all the neutrinos had exactly the same velocity vector, v u . Then 
Eq. (|) would give A ~ oj pl /^ 2 for the resonant modes and the dispersion relation 
(P complex solutions uj(k) with growth rates 7 ~ T 1 / 3 ^, around 10 10 s for the 
parameters shown above. This corresponds to the reactive instability put forward 
by Bingham et al. [|], ^ . 

However, that calculation is not realistic because the neutrinos do not have 
exactly the same direction of motion. As emphasized in ||, no matter how 
far the neutrinos are from the core, there is always an angular spread proportional 
to the neutrinosphere radius, R, 2a v ~ 2R/r. This causes a variation of u — k-v^ 
proportional to a u u) p i, orders of magnitude higher than any conceivable value for the 
width I7I of a resonance due to neutrino interactions. This means that u) p \ — k-v^ 
changes sign over the neutrino momentum distribution and only a vanishing small 
fraction of the neutrinos lie on the resonance. Such a situation is similar to Landau 
damping || |9| in an electron plasma and the imaginary part of A can be calculated 
by replacing (uj — k-v^) -1 with — m 8(uj p i — k-v^). It yields a neutrino contribution 
to 7 proportional to G 2 ? given by Eq. (|T5D, first obtained by Hardy and Melrose flC 



through other methods. The only possible exceptions, that we want to analyze here, 
are resonances located at some boundary of the phase space occupied by neutrinos. 
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One kinematic boundary is the velocity direction parallel to k (when k lies inside 
the neutrino velocity cone): it corresponds to the minimum angle between v„ and 
k and minimum value of u p i — kv„. The other boundary is the upper limit on 
the angle between w v and the radial direction, 9 < a v ~ R/r, due to the finite size 
of the neutrinosphere. Ideally, the distribution function would be discontinuous at 
9 = a u but that is not true as we will see. 

Consider first the case of a resonance at v v parallel to k i.e., io p \ = kv v . Then, 
the factor (1 — u; 2 /k 2 ) 2 in A vanishes unless the neutrinos are massive in which case 
v v ~ 1 — ml/2E%. Let v = 0J p i/k be the exact resonant speed. Then (1 — u; 2 /k 2 ) ~ 
ml/ El and u> p i — k-v is positive throughout the neutrino angular distribution. The 
problem is, the very mass that makes 1 — u; 2 /k 2 different from zero also causes a 
neutrino speed variation over the energy spectrum and a change of sign in u p i — k-v^. 
In fact, uo p i — kv v is negative for energies larger than Eq. Only energies obeying 
u p i — kv u < I7I can participate in the resonance of width |7|. The standard deviation 
AE U in the energy spectrum of Supernova neutrinos [11] , |T2|| is comparable to the 
average energy E v , AE V « E v /2.5, which implies a deviation Av v ~ 2m 2 /5.E 2 . In 
order that a significant fraction of the energy spectrum contributes to the resonance 
it is necessary that Av u < \j\/u! p i, but that puts an upper limit on the neutrino 
mass, ml /El < 5\j\/u p i, and, quite remarkably, on the factor 

ml < 5|t| 

That simply washes out the resonance because the numerator (1 — u; 2 /k 2 ) 2 in A 
becomes suppressed by 7 2 . To put it in another way, if there is a resonance for \ v 
parallel to k, the energy spectrum is always too broad to prevent uo p i — kv v from 
changing sign and departing from the resonance width. 

The other phase space boundary is the upper limit on the angle 9 between v„ and 
the radial direction, a v ~ R/r (R is the neutrinosphere radius and r the distance 
from the Supernova center). If the resonance lies on 9 = a u , i.e., uj p \ = k cos a u , then 
u p i — k v„ is essentially negative over the neutrino angular distribution. However, 
the distribution does not fall abruptly to zero at the polar angle a v . There are 
two reasons for this. First, the neutrinosphere radius depends on the v energy: 
the interaction cross sections increase with the v energy and consequently more 
energetic neutrinos suffer the last scattering in regions of lower density, farther from 
the center. Second, scattering is a statistical process by nature. Particles with 
the same energy suffer the last scattering at different radii according to a certain 
statistical distribution dependent upon the particular chemical and density profile 
of the medium (see Appendix). Both factors imply that the neutrinosphere has a 
considerable thickness and so has the neutrino angular aperture. This fact changes 
the way the distribution function depends on 9. 

Let Re and ARe be the average neutrinosphere radius and statistical uncertainty 
for neutrinos with well defined energy E v = E and = -Rg/r, = Ai?^/r the 



l-c 2 /k 2 ^<^. (6) 
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respective angular aperture and uncertainty. Assuming axial symmetry around the 
radial direction, the distribution function only depends on the energy and polar 
angle: f v = f v {E, 9). Its derivative with respect to 9 can be modeled as 

with g(E) ~ f u (E,0). It corresponds to a distribution function practically constant 
in the interval 9 < cue — 2 AaE and dropping to zero at 9 > ole + 2 AaE- In addition 
to AaE, the distribution function is also smoothed by the dependence of on the 
neutrino energy, which makes a total angular width Aa u = Aa E + daE/dE u AE U , 
centered on the polar angle a v . The consequence of this is that u p i — kv„ may change 
sign and depart from the resonance still within the angular boundary 9 pa a u ±2Aa u , 
if the resonance is not wide enough. 

Let vo be a particular vector of the resonance surface (defined by kv„ = uo p i) 
situated in the angular boundary i.e., with polar angle 9o close to a u , and azimutal 
angle O . The angular coordinates of a generic velocity vector v„ and wavevector k 
are denoted as (0,0) and (9 k ,(f) k ), respectively. Without loss of generality, (j) k = 0. 
The product k-Vj, and its variation from v are given by 

k-v^ = kv v (cos9 k cos 9 + sin9 k sin9 cos0) , (8) 
5k- v u = k e v 59 + k^v sin 9 5(f> , (9) 

where kg = k(— cos 9 k sin 9 + sin 9 k cos 9 Q cos O ) and k^ = —k sin 9 k sin O represent 
the components of k along the directions eg and e^, respectively. With k v = k-vo/^o 
they obey the identity kj+ + — 1. The angular displacement 

se = ± - 7C0S/j , (io) 

v/k% 2 - (k-v ) 2 

sin 9 50 = ± -. ^ sm ^ = (ii) 
v/k% 2 - (k-v ) 2 

with tan/? = k^/kg, causes a variation <5k-Vj, = ±7. Notice that k-v = uj p i and 
sin 0o c± sina^. If that displacement satisfies \59\ Aa v and |50| 7r, then 

— k -Vj, goes away to both sides of the resonance for velocity directions well 
inside the range 9 pa a u ± 2Aa u and therefore well inside the neutrino distribution. 
Of course, uj v \ — k-v F would have a definite sign if 9 > a v ± 2Aa v but then, the 
resonance would exist where the neutrino distribution function drops to zero i.e., 
there are no neutrinos. 

Our estimates of the neutrinosphere depth (see Appendix) give Aa v varying 
between 0.11 a u and 0.03 a u at different stages of Supernova evolution and a v = R/r 
between 0.1 and 0.03 at r = 300 km. On the other hand, I7I is orders of magnitude 
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below 10 9 u p i. This implies that the displacements calculated above satisfy 

\60\/Aa v < 10 3 | 7 |/KsinM , (12) 
\6<f>\/* < 10| 7 |/(^sin^) , (13) 

and therefore are well inside the neutrino angular boundary. The contrary would 
require a very small angle between k and v , Q^ v < 10 -6 for \~f\/ujpi < 10~ 9 , which 
leads to a huge suppression factor (1 — u 2 /k 2 ) 2 ~ Q\ v in A, which in turn further 
suppresses the values of 7 and so on, in other words, it falls in the case v parallel 
to k treated above in first place. 

The lesson from all this is that the contribution, j u , of neutrino weak interactions 
to the imaginary part of the wave frequency uj(k), is too small to produce resonances 
where the quantity A is inversely proportional to some power of 7„. On the contrary, 
A is independent of 7„ and 

lu = h u pl lm{A} . (14) 
Applying Landau prescription [|S], |Sj to Eq. (|5|), the result is 



uj p i 4a k 2 



(k 2 -u 2 pl ) 2 Jd 3 p u 5(u pl -k-v u )k-^ , (15) 



the same as obtained by Hardy and Melrose JED] from the study of stimulated emis- 
sion and absorption of plasmons by neutrinos. 7^ is proportional to Gp and sup- 
pressed by Gpn e n u /m e E u , down to ~ 10~ 28 even for electron and neutrino densities 
as high as n e = 10 29 cm -3 and n u ~ 10 30 cm -3 at 300 km from the Supernova center. 
The other point is, the rate 7^ does not drive the evolution of plasma waves because 
it is many orders of magnitude smaller than the damping rate caused by electron-ion 
collisions, j c , only two or three orders of magnitude below u p i for such high density 
plasmas 0, [|. 

More important than the growth rate is to know the energy transferred from 
neutrinos to the plasma. Keeping in mind that the overall wave growth is shut 
down by electron-ion collisional damping, the energy transferred per unity of time 
to a single mode is 2 , y u (k)ujpi(k). Assuming that the plasma waves obey a thermal 
equilibrium Bose-Einstein distribution, cut-off at the Debye wavenumber hp,, the 
total energy transferred per unity of time and volume is 

• 9 f d3k lu(k)uj p i 

P = 2 J (2,)3 e ^-r (16) 

Integrating in time and volume and dividing by the total neutrino energy, 47rr 2 n u E u At, 
one obtains the fraction of total energy transferred per neutrino, AE / E u = fdrp/ E v n v 

For the sake of argument let us assume that some new interaction produces 
a contribution like A in Eq. (||), but without the factor (1 — cj^/k 2 ) 2 , so that a 
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resonance is possible for k parallel to \ v {u p i = k) making A oc 7 _1 . The collisionless 
dispersion relation ([I]) gives a growth rate 7 oc T 1 I 2 uj p i that is proportional to G-p 
rather than G F . However, these resonant modes are limited to the very thin shell 
k = ujpi ±7 in k space. That makes an extra factor of 7 in the integration of Eq. flTED 
and the transferred energy goes as 7 2 , proportional to G F not to Gp. Taking only 
the main factors, one obtains AE/E U ~ G\T e io^ vl E~ 2 r , about 10~ 18 for the same 
numbers used before and T e = 100 KeV. This is very far from the few percent 
needed for a Supernova explosion. It discourages the application of non-standard 
weak interactions seeking for non-linear resonant effects. 



3 Conclusions 

We analyzed different instabilities that could possibly emerge from the interaction 
between neutrinos and electron plasma waves in a Supernova environment. The hy- 
pothetical neutrino induced resonances are too narrow to embrace the total angular 
spread of the neutrino stream due to the finite size of the neutrinosphere 0, 
But they are also too narrow to contain the boundary regions of the neutrino phase 
space distribution. When the resonant neutrino velocity vector is parallel to the 
wavevector, a neutrino mass is needed and the energy spectrum is too broad to keep 
the neutrino speed, v u ~ 1 — mr v /2E 2 U1 inside the resonance. When the resonant 
velocities are in the boundary of the velocity angular distribution, the depth of the 
neutrinosphere is too large to make an angular boundary abrupt enough. In both 
cases the boundaries of the neutrino phase space distribution are too wide to prevent 
a departure from the poles to both sides of the resonances. As a result, non-linear 
resonant effects are not possible and the growth or damping rates are always linear 
in the neutrino flux and G F . They correspond to a balance between stimulated 



Cerenkov emission and absorption of plasma waves by neutrinos ||10|| , analogous to 
Landau damping. The energy that could possibly be transferred from neutrinos to 
plasma waves seems to be vanishing small even for non-standard neutrino interac- 
tions. 
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Appendix: Neutrinosphere depth 



Neutrinos like any other form of radiation are not emitted from an ideally surface but 
from a shell with a certain depth. The neutrino 'optical' depth at a point of radius 
r can be defined as the Roseland depth, r(r) = J°° dr/X(r), where A -1 = kp = an 
is the inverse mean free path, a function of the opacity k and density p, or cross 
section o and number density of target particles, n. r(r) is a measurement of the 
number of collisions the neutrinos suffer moving from the radius r to infinity. The 
neutrinosphere can be arbitrarily defined as the surface, of radius R, where r = 2/3, 
but the bulk of neutrino last scatterings spreads between the r = 1 and r = 1/3 
surfaces. The way this translates as a radial distribution depends on the chemical 
and density profile of the medium. For an exponentially decreasing density near the 
neutrinosphere, p(r) oc lCT r//£ , the radial separation between r = 1 and r = 1/3 is 
approximately 2AR ml/2. Numerical simulations [13] show that the length scale 
I varies from around R/5 to R/20 at different stages of Supernova. Then, AR is 
about R/20 to R/80. 

The other source of radial spread is the opacity or cross section dependence 
on the neutrino energy [12|], typically k oc E%. At constant 'optical' depth r, 
the neutrinosphere radius increases with the energy at a rate that we estimate as 
5R/5E y = 21/ E v In 10. Since the neutrino energy spectrum [TIL LJ3] has a standard 



deviation AE U m E v /2.h we obtain a radial deviation AR m 1/3 w R/15 — R/60. 
The joint effect of statistical fluctuations and energy spectrum is a neutrinosphere 
depth varying between 2AR m 0.23R and 0.06-R. This implies a finite width in the 
boundary of the neutrino angular aperture at large distances from the Supernova 
core as discussed in the text. 
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